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Abstract 

We conjecture that because of color confinement, the physical vacuum forms an event 
horizon for quarks and gluons which can be crossed only by quantum tunneling, i.e., 
through the QCD counterpart of Hawking radiation by black holes. Since such radiation 
cannot transmit information to the outside, it must be thermal, of a temperature de- 
termined by the chromodynamic force at the confinement surface, and it must maintain 
color neutrality. We explore the possibility that the resulting process provides a com- 
mon mechanism for thermal hadron production in high energy interactions, from e + e~ 
annihilation to heavy ion collisions. 
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1 Introduction 



The aim of this paper is to develop a conceptual framework for a universal form of thermal 
multihadron production in high energy collisions. Our work is based on two seemingly 
disjoint observations: 

• Color confinement in QCD does not allow colored constituents to exist in the physi- 
cal vacuum, and thus in some sense creates a situation similar to the gravitational 
confinement provided by black holes. 

• Numerous high energy collision experiments have provided strong evidence for the 
thermal nature of multihadron production, indicating a universal hadronization tem- 
perature T H ~ 150 - 200 MeV. 

We want to suggest that quantum tunnelling through a color event horizon, as QCD 
counterpart of Hawking-Unruh radiation from black holes, can relate these observations 
in a quite natural way. 

The idea that the color confinement of quarks and gluons in hadrons may have a dual 
description in terms of a theory in curved space-time is not new. Both gravitational 
confinement of matter inside a black hole [1] and the de Sitter solution of the Einstein 
equations with a cosmological constant describing a "closed" universe of constant cur- 
vature [2] have been proposed as possible descriptions of quark confinement. Soon it 
became clear that asymptotic freedom [3] and the scale anomaly [4,5] in QCD completely 
determine the structure of low-energy gluodynamics [6]. This effective theory can be con- 
veniently formulated in terms of the Einstein-Hilbert action in a curved background. At 
shorter distances (inside hadrons), the effective action has the form of classical Yang-Mills 
theory in a curved (but conformally flat) metric [7]. The "cosmological constant" present 
in this theory corresponds to the non-perturbative energy density of the vacuum, or the 
"gluon condensate" [8]. 

It is worthwhile to also mention here the well-known conjectured holographic correspon- 
dence between the large N limit of supersymmetric Yang-Mills theory in 3 + 1 dimensions 
and supergravity in an anti-de Sitter space-time sphere, AdS 5 x S 5 [9]. This example 
illustrates the possible deep relation between Yang-Mills theories and gravity; however, 
a conformal theory clearly differs from the examples noted above, in which the scale 
anomaly (describing the breaking of conformal invariance by quantum effects) was used 
as a guiding principle for constructing an effective curved space-time description. 

Let us assume that color confinement indeed allows a dual description in terms of the 
gravitational confinement of matter inside black holes. What are the implications of 
this hypothesis for hadronic physics? Hawking [10] showed that black holes emit thermal 
radiation due to quantum tunneling through the event horizon. Shortly afterwards, Unruh 
[11] demonstrated that the presence of an event horizon in accelerating frames also leads 
to thermal radiation. It was soon conjectured that the periodic motion of quarks in a 
confining potential [12], or the acceleration which accompanies inelastic hadronic collisions 
[13-15], are associated with an effective temperature for hadron emission. 

Recently, a QCD-based picture of thermal production based on the parton description of 
high energy hadronic collisions has been proposed [16, 17]. The effective temperature T 
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is in this case determined either by the string tension a, with a relation 




(1) 



or, in the gluon saturation regime, by the saturation momentum Q s describing the strength 
of the color fields in the colliding hadrons or nuclei, with T ~ (Q s /2n). 

Turning now to the second observation, we recall that over the years, hadron production 
studies in a variety of high energy collision experiments have shown a remarkably universal 
feature. From e + e~ annihilation to p — p and p — p interactions and further to collisions 
of heavy nuclei, with energies from a few GeV up to the TeV range, the production 
pattern always shows striking thermal aspects, connected to an apparently quite universal 
temperature around T H ~ 150 — 200 MeV [18]. As a specific illustration we recall that 
the relative abundance of two hadron species a and b, of masses m a and m&, respectively, 
is essentially determined by the ratio of their Boltzmann factors [19], 



What is the origin of this thermal behaviour? While high energy heavy ion collisions 
involve large numbers of incident partons and thus could allow invoking some "thermal- 
isation" scheme through rescattering, in e + e~ annihilation the predominant initial state 
is one energetic qq pair, and the number of hadronic secondaries per unit rapidity is too 
small to consider statistical averages. The case in p — p/p — p collisions is similar. 

This enigma has led to the idea that all such collision experiments result in the formation 
of a strong color field "disturbing" the physical vacuum. The disturbed vacuum then 
recovers by producing hadrons according to a maximum entropy principle: the actually 
observed final state is that with the largest phase space volume. While this provides an 
intuitive basis for a statistical description, it does not account for a universal temperature. 
Why don't more energetic collisions result in a higher hadronization temperature? 

A further piece in this puzzle is the observation that the value of the temperature deter- 
mined in the mentioned collision studies is quite similar to the confinement /deconfinement 
transition temperature found in lattice studies of strong interaction thermodynamics [20]. 
While hadronization in high energy collisions deals with a dynamical situation, the ener- 
gy loss of fast color charges "traversing" the physical vacuum, lattice QCD addresses 
the equilibrium thermodynamics of unbound vs. bound color charges. Why should the 
resulting critical temperatures be similar or even identical? 

We shall here consider these phenomena as reflections of the QCD counterpart of the 
Hawking radiation emitted by black holes [10]. These ultimate stellar states provide a 
gravitational form of confinement and hence, as already noted, their physics was quite 
soon compared to that of color confinement in QCD [1,2], where colored constituents 
are confined to "white holes" (colorless from the outside, but colored inside). It should 
be emphasized from the outset that in contrast to the original black hole physics in 
gravitation, where confinement is on a classical level complete, in QCD confinement refers 
only to color-carrying constituents; thus, e.g., photons or leptons are not affected. 

In black hole physics, as noted above, it was shown that the event horizon for systems 
undergoing uniform acceleration leads to quantum tunnelling and hence to thermal radia- 
tion [11]. Our aim here is to show that such Hawking-Unruh radiation, as obtained in the 



R(a/b) ~ exp -{(m a - m h )/T H }. 



(2) 
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specific situation of QCD, provides a viable account for the thermal behavior observed 
in multihadron production by high energy collisions. Furthermore, in the process we also 
want to elucidate a bit the common origin of the "limiting temperature" concepts which 
have arisen in strong interaction physics over the years. 

We begin by reviewing those features of black hole physics and Hawking radiation which 
are relevant for our considerations, and then discuss how they can be implemented in 
QCD. In particular, we show that modifications of the effective space-time structure, in a 
perturbative approach as well as in a non-perturbative treatment based on a large-scale 
dilaton field, lead to an event horizon in QCD. 

Following this, we present the main conceptual consequences of our conjecture. 

• Color confinement and the instability of the physical vacuum under pair produc- 
tion form an event horizon for quarks, allowing a transition only through quantum 
tunnelling; this leads to thermal radiation of a temperature Tq determined by the 
string tension. 

• Hadron production in high energy collisions occurs through a succession of such 
tunnelling processes. The resulting cascade is a realization of the same partition 
process which leads to a limiting temperature in the statistical bootstrap and dual 
resonance models. 

• The temperature Tq of QCD Hawking-Unruh radiation can depend only on the 
baryon number and the angular momentum of the deconfined system. The former 
could provide a dependence of Tq on the baryon number density, while the angular 
momentum pattern of the radiation allows a centrality-dependence of Tq and elliptic 
flow. 

• In kinetic thermalization, the initial state information is successively lost through 
collisions, converging to a time-independent equilibrium state. In contrast, the 
stochastic QCD Hawking radiation is "born in equilibrium", since quantum tun- 
nelling a priori does not allow information transfer. 

2 Event Horizons in Gravitation and in QCD 
2.1 Black holes 

A black hole is formed as the final stage of a neutron star after gravitational collapse [24]. 
It has a mass M concentrated in such a small volume that the resulting gravitational field 
confines all matter and even photons to remain inside the event horizon R of the system: 
no causal connection with the outside is possible. As a consequence, black holes have 
three (and only three) observable properties: mass M, charge Q and angular momentum 
J. This section will address mainly black holes with Q = J = 0; we shall come back to 
the more general properties in section 4. We use units of h — c — 1. 

The event horizon appears in a study of the gravitational metric, which in flat space has 
the form 

ds 2 = gdt 2 - g' 1 dr 2 - r 2 [d6 2 + sin 2 6d<f) 2 } , (3) 
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using units where c = 1. The field strength of the interaction is contained in the coefficient 

g(r), 

. ( 2GM\ 
9(r) = (1- — ) , (4) 

leading back to the Minkowski metric in the large distance limit r — > oo. The vanishing 
of g(r) specifies the Schwarzschild radius R as event horizon, 

R = 2GM. (5) 

It is interesting to note that the mass of a black hole thus grows linearly with R, analogous 
to the behaviour of the confining potential in strong interactions: M(R) = (2G) _1 R. 

Classically, a black hole would persist forever and remain forever invisible. On a quan- 
tum level, however, its constituents (photons, leptons and hadrons) have a non-vanishing 
chance to escape by tunnelling through the barrier presented by the event horizon. Equiv- 
alently, we can say that the strong force field at the surface of the black hole can bring 
vacuum fluctuations on-shell. The resulting Hawking radiation [10] cannot convey any 
information about the internal state of the black hole; it must be therefore be thermal. 
For a non-rotating black hole of vanishing charge (denoted as Schwarzschild black hole), 
the first law of thermodynamics, 

dM = TdS (6) 
combined with the area law for the black hole entropy [25], 

S = 7T < 7 > 

leads to the corresponding radiation temperature 

This temperature is inversely proportional to the mass of the black hole, and since the 
radiation reduces the mass, the radiation temperature will increase with time, as the black 
hole evaporates. For black holes of stellar size, however, one finds T B h < 2x 1CT 8 °K, 
which is many orders of magnitude below the 2.7 °K cosmic microwave background, and 
hence not detectable. 

It is instructive to consider the Schwarzschild radius of a typical hadron, assuming a mass 
m ~ 1 GeV: 

R had _ 1 3 x 1Q -38 Q e y-i ~ 2.7 x 10~ 39 fm. (9) 

To become a gravitational black hole, the mass of the hadron would thus have to be 
compressed into a volume more than 10 100 times smaller than its actual volume, with a 
radius of about 1 fm. On the other hand, if instead we increase the interaction strength 
from gravitation to strong interaction [1], we gain in the resulting "strong" Schwarzschild 
radius Rg ad a factor 
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where a s is the dimensionless strong coupling and Gm 2 the corresponding dimensionless 
gravitational coupling for the case in question. This leads to 

s ~ m ' K ' 

with the effective value of a s ~ 0(1) we thus get R^ ad ~ 0(1) fmEl In other words, 
the confinement radius of a hadron is about the size of its "strong" Schwarzschild radius, 
so that we could consider quark confinement as the strong interaction version of the 
gravitational confinement in black holes [1,2]. 

We had seen that the mass of a black hole grows linearly with the event horizon, M = 
(1/2G)R, so that in gravitation 1/2G plays the role of the string tension in strong inter- 
action physics. The replacement GM 2 — > a s here leads to 

o 

777 

g — ~ 0.16 GeV 2 , (12) 

2a s 



<s 



if one uses the mentioned effective saturation value a s ~ 3 [26] . The value of a s ~ 1 thus 
gives a reasonable string tension as well as a reasonable radius. 



2.2 Quasi- Abelian case 

The appearence of an event horizon occurs in general relativity through the modification 
of the underlying space-time structure by the gravitational interaction. Such modifica- 
tions have also been discussed for other interactions. In particular, it was noted that 
in electrodynamics, non-linear in-medium effects can lead to photons propagating along 
geodesies which are not null in Minkowski space-time; this can even lead to photon trap- 
ping, restricting the motion of photons to a compact region of space [27]. Thus, an 
effective Lagrangian C(F) depending on a one-parameter background field, F = F^F^", 
results in a modified metric 

9^u = ri^C - 4F a ^F«£", (13) 

where the primes indicate first and second derivatives with respect to F. Hence 

g 00 =£'- AFC" = (14) 

defines the radius of the compact region of the theory, i.e., the counterpart of a black 
hole [27]. 

QCD is an inherently non-linear theory, with the physical vacuum playing the role of a 
medium [28]. The general structure of the effective Lagrangian in a background field F, 
compatible with gauge invariance, renormalization group results [3] and trace anomaly, 
has the unique form [29] 

Cqcd = \ F ^j^ = \V V <9F). (15) 

*In fact, some studies [26] indicate that at large distances, the strong coupling freezes at a s ~ 3; in 
that case the corresponding radius becomes R^ ad ~ 1 fm. 
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Here e(gF) is the dielectric "constant" of the system in the presence of the background 
field; the F-dependence of e(F) effectively turns the QCD vacuum into a non-linear 
medium. On a one-loop pertubative level we have 

£ ( 9 F)-l-A>(Qln^, (16) 

where (3q = (11N C — 2Nf)/A8ir 2 , with N c and Nf specifying the number of colors and 
flavors, respectively. Using this form in the formalism of [27] leads to g t changing sign 
(i.e., g t = 0) at 

gF* = A 2 ex P {-4n/(3 g 2 }, (17) 

indicating a possible horizon at r* ~ 1 / y/gF* [30] . It is clear that this line of argument 
can at best provide some hints, since we used the lowest order perturbative form of the 
beta-function, even though at the horizon perturbation theory will presumably break 
down. Nevertheless, we believe that it suggests the possibility of an event horizon for 
QCD; the crucial feature is the asymptotic freedom of QCD [3], which leads to e < 1 and 
allows goo to vanish even without the external medium effects required in QED. 



2.3 Non-Abelian case 



Indeed, a different and more solid suggestion that in QCD there is an event horizon comes 
from studying the theory on a curved background. For gluodynamics, such a program is 
discussed in [7]. Classical gluodynamics is a scale- invariant theory, but quantum fluctu- 
ations break this invariance, with the trace of the energy-momentum tensor introducing 
non-perturbative effects, associated with the vacuum energy density e vac . It was shown [6] 
that low-energy theorems can be used to determine the form of the effective Yang-Mills 
Lagrangian in a curved but conformally flat metric 



x ) 



'It 



where the dilaton field h(x) is coupled to the trace of the energy momentum tensor, 9£. 
The resulting action has the form 



S 



d x 



4e 



1 



l^e\d»hY-UF« v y + e 2h {e™ 



3 m 



G 



4 



7 fj,\pert)- 



(19) 



here e vac is the absolute value of the energy density of the vacuum and mg the dilaton 
mass; the trace of the energy-momentum tensor has been separated into perturbative and 
non-perturbative contributions, 
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9£\pert+ < 8% >= ^ | pert 



4e^ 



(20) 



The crucial point for our considerations is that the first term of eq. (Tl9~|) can be written 

as 

Ry/=9, (21) 



3 

defining R as the Ricci scalar of the theory. Hence eq. (JT9l has the structure of an 
Einstein-Hilbert Lagrangian of gluodynamics in the presence of an effective gravitation, 



G 



d x 



a .2 _|_ g2/i ^g vac 



1 li\ pert)- 



(22) 
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where G is now given by 

1 64vre vac 



G 3 m 2 G 



(23) 



The relation 1/2G — > a between G and the string tension conjectured above then leads 
to 



vac 



32vr e 

° = -5 — r- 24 

3 7T7.Q 

On the other hand, the string tension is just the energy density of the vacuum times the 
transverse string area, 

a = e vac 7rr 2 . (25) 
Combining relations ff24l) and (|25p . we have 

/S2~ 1 

r T = W ~ 0.4 fin, (26) 

V 3 m G 

using me — 1.5 GeV for the scalar glueball mass. Eq. (126]) thus gives us the transverse 
extension or horizon of the string. 

From eqs. ( l24l) or (1251) we can obtain a further consistency check. Given the glueball mass 
and the string tension a ~ 0.16 GeV 2 , we find for the vacuum energy density 

e vac _ JL CT m 2 Q ^ Q.013 GeV 4 ~ 1.7 GeV/fm 3 . (27) 

327T 

This is the value for pure gluodynamics; since the energy density is related to the trace 
of the energy-momentum tensor by the relation fTSOj) . and 



M ( F - )i~- b ^ 
2g y " u) 32tt 2 K ^ vl ' 



K = ^ ( Fa J * {FL) 2 , (28) 



with the coefficient b = HN c — 2Nf of the /^-function, we can estimate that for three-flavor 
QCD 

ll N c -2N f cVac = 9_ 

UNn ' 11" 



/V — 2 \ t L ) 
e ^ D = ±— — l e vac = — e vac ~ 0.01 GeV 4 , (29) 



which is in perfect agreement with the original value of the gluon condensate [8] 

— G 2 ) ~ 0.012 GeV 4 ; (30) 

7T / 

note that e^g D = 27/32 ((a s /n)G 2 )). 



3 Hyperbolic Motion and Hawking- Unruh Radiation 

In general relativity, the event horizon appeared as consequence of the geometrized grav- 
itational force, but its occurrence and its role for thermal radiation was soon generalized 
by Unruh [11]. A system undergoing uniform acceleration a relative to a stationary ob- 
server eventually reaches a classical turning point and thus encounters an event horizon. 
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Let us recall the resulting hyperbolic motion [34]. A point mass m subject to a constant 
force F satisfies the equation of motion 



d 



dty/l 



mv 



(31) 



where v(t) = dx/dt is the velocity, normalized to the speed of light c = 1. This equation 
is solved by the parametric form through the so-called Rindler coordinates, 



x 



£ cosh ar t = £ sinh ar, 



(32) 



where a = F/m denotes the acceleration in the instantaneous rest frame of m, and r 
the proper time, with dr = yl — v 2 dt. If we impose the boundary condition that the 
velocity at t = vanishes, we have £ = 1/a and x(t = 0) — 1/a. The resulting world line 
is shown in Fig. [TJ It corresponds to the mass m coming from x = oo at t = — oo at with 
a velocity arbitrarily close to that of light, decelerating uniformly until it comes to rest 
at the classsical turning point xh = —(l/a),t = 0. Subsequently, it accelerates again and 
returns to x = oo at t = oo, approaching the speed of light. For given a, the light cone 
originating at a distance xh = 1/a away from the turning point of m defines a space-time 
region inaccessible to m: no photon in this region can (classically) ever reach m, in much 
the same way as photons cannot escape from a black hole. Here the acceleration is crucial, 
of course; if m stops accelerating, it will eventually become visible in the "hidden region". 



event horizon 



hiddei 
regioi 


1 

1 // 

1 mass m 







Figure 1: Hyberbolic motion 



The metric of such a accelerating system becomes in spherical coordinates [31] 

ds 2 = i 2 a 2 dr 2 - d£ 2 - £ 2 cosh 2 ar(d0 2 + sin 2 6d<f) 2 ), (33) 



(34) 



which we want to compare to the black hole metric ([3]). Making in the latter the coordinate 
transformation [32] 

V = 

K 

where the surface gravity k is given by 



K 




(35) 



r=R 



we obtain for r 



R the black hole form 

ds 2 = V 2 n 2 dt 2 - d V 2 - R 2 (d6 2 + sin 2 



(36) 



When we compare eqs. (1331) and (136]) . it is evident that the system in uniform acceleration 
can be mapped onto a spherical black hole, and vice versa, provided we identify the surface 
gravity k with the acceleration a. 

The vacuum through which m travels is, for a stationary observer, empty space. On a 
quantum level, however, it contains vacuum fluctuations. The accelerating mass m can 
bring these on-shell, using up a (small) part of its energy, so that for m the vacuum 
becomes a thermal medium of temperature 



T, 



u 



a 

2tt" 



(37) 



Consider such a fluctuation into an e + e~ pair, flying apart in opposite directions. One 
electron is absorbed by the mass m, the other penetrates into the "hidden region" and can 
never be detected by m (see Fig. [2j Since thus neither an observer on m nor a stationary 
observer in the hidden region can ever obtain access to full information, each will register 
the observed radiation as thermal (Einstein- Podolsky- Rosen effect [33]). In other words: 
the accelerating mass m sees the vacuum as a physical medium of temperature Tu, while 
a stationary observer in the hidden region observes thermal radiation of temperature Tjj 
as a consequence of the passing of m. 



e 


/ mass m 







Figure 2: Unruh radiation 



We here also mention that the entropy in the case of an accelerating system again becomes 
1/4 of the event horizon area, as in the black hole case, so that also here the correspondence 
remains valid [35]. 

In the case of gravity, we have the force 



ma 



G 



mM 



(38) 



on a probe of mass m. With R = 2 GM for the (Schwarzschild) black hole radius, we have 
a = 1/(4 GM) for the acceleration at the event horizon and hence the Unruh temperature 
(1371) leads back to the Hawking temperature of eq. (JHJ). 
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In summary, we note that constant acceleration leads to an event horizon, which can be 
surpassed only by quantum tunnelling and at the expense of complete information loss, 
leading to thermal radiation as the only resulting signal. 



4 Pair Production and String Breaking 

In the previous section, we had considered a classical object, the mass m, undergoing 
accelerated motion in the physical vacuum; because of quantum fluctuations, this vacuum 
appears to m as a thermal medium of temperature Ty. In this section, we shall first 
address the modifications which arise if the object undergoing accelerated motion is itself 
a quantum system, so that in the presence of a strong field it becomes unstable under 
pair production. Next we turn to the specific additional features which come in when 
the basic constituents are subject to color confinement and can only exist in color neutral 
bound states. 

As starting point, we consider two-jet e + e~ annihilation at cms energy y/s, 

e + e~ — > 7* — > qq — > hadrons. (39) 

The initially produced qq pair flies apart, subject to the constant confining force given 
by the string tension a; this results in hyperbolic motion [14] of the type discussed in the 
previous section. At t = 0, the q and q separate with an initial velocity Vq = pj \Jp 2 + m 2 , 
where p ~ \fsj2 is the momentum of the primary constituents in the overall cms and m 
the effective quark mass. We now have to solve Eq. (13*1]) with this situation as boundary 
condition; the force 

F = a, (40) 
is given by the string tension a binding the qq system. The solution is 



x = [1 - VI -vot + t*] (41) 
with x = x/xq and t = t/xo; here the scale factor 

mil , . 

x = == = - 7 42 

° Jl^l a 



is the inverse of the acceleration a measured in the overall cms. The velocity becomes 



„„ = ffe ("°/2)-f . (43) 
dt J l- v t + t 2 



it vanishes for 
thus defining 

<7 \ ) 2a 

as classical turning point and hence as the classical event horizon measured in the overall 
cms (see Fig. [3]) . 



f. = |^ t . = |^ 7 , (44) 
z la 

x(f ) = - 7 (l - Jl - (v 2 /A)) ~ # (45) 
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Xq X(t*) X 



Figure 3: Classical and quantum horizons in qq separation 

Eq. (|45p allows the q and the q to separate arbitrarily far, provided the pair has enough 
initial energy; this clearly violates color confinement. Our mistake was to consider the 
qq system as classical; in quantum field theory, it is not possible to increase the potential 
energy of a given qq state beyond the threshold value necessary to bring a virtual qq pair 
on-shell. In QED, the corresponding phenomenon was addressed by Schwinger [36], who 
showed that in the presence of a constant electric field of strength 8 the probability of 
producing an electron-positron pair is given by 

P(M, 8) ~ exp{-nm 2 Je8 }, (46) 

with m e denoting the electron mass and e denoting the electric charge. This result is 
in fact a specific case of the Hawking-Unruh phenomenon, as shown in [16]. In QCD, 
we expect a similar effect when the string tension exceeds the pair production limit, i.e., 
when 

a x >2m (47) 

where m specifies the effective quark mass. Beyond this point, any further stretching of 
the string is expected to produce a qq pair with the probability 

P(M, a) ~ exp{-nm 2 /cx}, (48) 

with the string tension a replacing the electric field strength e8. This string breaking 
acts like a quantum event horizon x q = 2 m/ a, which becomes operative long before 
the classical turning point is ever reached (see Fig. [3]). Moreover, the resulting allowed 
separation distance for our qq pair, the color confinement radius xq, does not depend on 
the initial energy of the primary quarks. 

There are some important differences between QCD and QED. In case of the latter, the 
initial electric charges which lead to the field 8 can exist independently in the physical 
vacuum, and the produced pair can be simply ionized into an e + and an e~ . In contrast, 
neither the primary quark nor the constituents of the qq pair have an independent ex- 
istence, so that in string breaking color neutrality must be preserved. As a result, the 
Hawking radiation in QCD must consist of qq pairs, and these can be produced in an in- 
finite number of different excitation states of increasing mass and degeneracy. Moreover, 
the qq pair spectrum is itself determined by the strength a of the field, in contrast to the 
exponent m 2 e /8 in eq. (JIB"]) , where the value of 8 has no relation to the electron mass m e . 
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Hadron production in e + e annihilation is believed to proceed in the form of a self-similar 
cascade [37,38]. Initially, we have the separating primary qq pair, 

7 [qq] (49) 

where the square brackets indicate color neutrality. When the energy of the resulting 
color flux tube becomes large enough, a further pair q\q\ is excited from the vacuum by 
two-gluon exchange (see Fig. 0J , 

7 -> \q[Mi]q}- (50) 

Although the new pair is at rest in the overall cms, each of its constituents has a transverse 
momentum determined, through the uncertainty relation, by the transverse dimension 
rrp of the flux tube. The slow q\ now screens the fast primary q from its original partner 
q, with an analoguous effect for the gi and the primary antiquark. To estimate the qq 
separation distance at the point of pair production, we recall that the thickness of the 
flux tube connecting the qq pair is in string theory given by [41] 

2 K 1 

r 2 = jf_\-_J_ (51) 

k=0 

where K is the string length in units of an intrinsic vibration measure. Lattice studies [43] 
show that for strings in the range of 1 - 2 fm, the first string excitation dominates, so that 
we have 

/ 2 

r T = c \ — , (52) 

V 7TCT" 

with Co — 1 or slightly larger. Higher excitations lead to a greater thickness and eventually 
to a divergence (the "roughening" transition). From the uncertainty relation we then have 

kr = • (53) 

C V z 

With this transverse energy is included in eq. (JSj), we obtain for the pair production 
separation xq 



ax q = 2Jm 2 + k? =J> x t r ^ 



<-i 



-Jw? + (7ra/2c 2 ) ~ ~ 1 fm, (54) 

o y ctCq 



with cr = 0.2 GeV , m <C ex, and c ~ 1 



Once the new pair is present, we have a color-neutral system qqiqiq; but since there is a 
sequence of connecting string potentials qq~i, q±qi and gig, the primary string is not yet 
broken. To achieve that, the binding of the new pair has to be overcome, i.e., the gi has to 
tunnel through the barrier of the confining potential provided by q~i, and vice versa. Now 
the q excerts a longitudinal force on the g~i, the g on the gi, resulting in a longitudinal 
acceleration and ordering of gi and q\. When (see Fig. H|) 



ax(q 1 q 1 ) = 2Jm 2 + k\, (55) 
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Figure 4: String breaking through qq pair production 

the qi reaches its qiq~i horizon; on the other hand, when 

ax(qqi) = 2\Jm? + kj-, (56) 

the new flux tube qq\ reaches the energy needed to produce a further pair g 2 g 2 . The 
g 2 screens the primary q from the q± and forms a new flux tube gg 2 . At this point, the 
original string is broken, and the remaining pair gig 2 form a color neutral bound state 
which is emitted as Hawking radiation in the form of hadrons, with the relative weights 
of the different states governed by the corresponding Unruh temperature. The resulting 
pattern is schematically illustrated in Fig. HI 

To determine the temperature of the hadronic Hawking radiation, we return to the original 
pair excitation process. To produce a quark of momentum kx, we have to bring it on-shell 
and change its velocity from zero to v = kx/(m 2 + A;^) 1 / 2 ~ 1. This has to be achieved in 
the time of the fluctuation determined by the virtuality of the pair, At =1/ AE ~ \/2kr- 
The resulting acceleration thus becomes 

a = — ^ ~ 2 k T ~ V2^a 7c ~ 1 GeV, (57) 
At 

which leads to 

T Q = — ~ — \f~~ ^ 160 - 180 MeV (58) 
2ir c V 2tt 

for the hadronic Unruh temperature. It governs the momentum distribution and the 
relative species abundances of the emitted hadrons. 

A given step in the evolution of the hadronization cascade of a primary quark or antiquark 
produced in e + e~ annihilation thus involves several distinct phenomena. The color field 
created by the separating q and q produces a further pair q x q\ and then provides an 
acceleration of the qi, increasing its longitudinal momentum. When it reaches the q\q~i 
confinement horizon, still another pair g 2 ?2 is excited; the state gig 2 is emitted as a 
hadron, the g 2 forms together with the primary q a new flux tube. This pattern thus 
step by step increases the longitudinal momentum of the "accompanying" as well as 
of the emitted hadron. This, together with the energy of the produced pairs, causes a 
corresponding deceleration of the primary quarks q and q, in order to maintain overall 
energy conservation. In Fig. [5j we show the world lines given by the acceleration qi 
(qi — > qi+i) and that of formation threshold of the hadrons qiqi+i and the corresponding 
antiparticles. 
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Figure 5: Quark acceleration and hadronization world lines 



The energy loss and deceleration of the primary quark q in this self-similar cascade, 
together with the acceleration of the accompanying partner from the successive pairs 
brings q and closer and closer to each other in momentum, from an initial separation 
qq~i of y/s/2, until they finally are combined into a hadron and the cascade is ended. The 
resulting pattern is shown in Fig. [6j 




The number of emitted hadrons, the multiplicity v(s), follows quite naturally from the 
picture presented here. The classical string length, in the absence of quantum pair for- 
mation, is given by the classical turning point determined in eq. fj45|) . The thickness of a 
flux tube of such an "overstretched" string is known [41]; from eq. (jSTl) we get 

R 2 t = — E^T^— 1H2JT, (59) 

1XO 2k + 1 7TCT 

where K is the string length. From eq. (1451) we thus get 

R 2 T ~ — In (60) 

7TO" 

for the flux tube thickness in the case of the classical string length. In parton language, the 
logarithmic growth of the transverse hadron size is due to parton random walk (" Gribov 
diffusion" [42]); this phenomenon is responsible for diffraction cone shrinkage in high- 
energy hadron scattering. 

Because of pair production, the string breaks whenever it is stretched to the length x q 
given in eq. ( l54|) ; its thickness tt at this point is given by eq. (|5Tjl . The multiplicity can 
thus be estimated by the ratio of the corresponding classical to quantum transverse flux 
tube areas, 

R 2 

u(s) ~ — 5- ~ In y/s, (61) 
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and is found to grow logarithmically with the e + e annihilation energy, as is observed 
experimentally over a considerable range. 

We note here that in our argumentation we have neglected parton evolution, which would 
cause the emitted radiation (e.g., qiq2 in Fig. HJ) to start another cascade of the same 
type. Such evolution effects result eventually in a stronger increase of the multiplicity. 
The formation of a white hole does not affect the production of hard processes at early 
times (e.g., multiple jet production), which is responsible for an additional growth of the 
measured multiplicity. 

A further effect we have not taken into account here is parton saturation. At sufficiently 
high energy, stronger color fields can lead to gluon saturation and thus to a higher temper- 
ature determined by the saturation momentum [16]. The resulting system then expands 
and hadronizes at the universal temperature determined by the string tension. 

It interesting to compare the separation of two energetic light quarks, as we have consid- 
ered here, with that of two static heavy quarks Q and Q. From quarkonium studies it is 
known that 

2(M D - m c ) ~ 2{M B - m b ) ~ 1.2 GeV, (62) 

where Mjj(Mb) and m c {mh) are the masses of open charm (beauty) mesons and of the 
corresponding charm (beauty) quarks, respectively. The energy needed to separate a 
heavy QQ pair thus is independent of the mass of the heavy quarks, indicating that 
the string breaking involved here is really a consequence of the vacuum, through qq pair 
excitation. With 

ax Q ~ 1.2 GeV x Q ~ 1.2 fm (63) 

we find that the resulting separation threshold for pair excitation agrees well with that 
found above in eq. (J51j) . Lattice QCD studies lead to similar results. 

Up to now, we have considered hadron production in e + e~ annihilation, in which the 
virtual photon produces a confined colored qq pair as a "white hole". Turning now to 
hadron-hadron collisions, we note that here two incident white holes combine to form a 
new system of the same kind, as schematically illustrated in Fig. [7J Again the resulting 
string or strong color field produces a sequence of qq pairs of increasing cms momentum, 
leading to the well-known multiperipheral hadropro duct ion cascade shown in Fig. [SJ We 
recall here the comments made above concerning parton evolution and saturation; in 
hadronic collisions as well, these phenomena will affect the multiplicity, but not the relative 
abundances. 

In the case of heavy ion collisions, two new elements enter. The resulting systems could 
now have an overall baryon number, up to B = 400 or more. To take that into account, 
we need to consider the counterpart of charged black holes. Furthermore, in heavy ion 
collisions the resulting hadron production can be studied as function of centrality, and 
peripheral collisions could lead to an interaction region with an effective overall angular 
momentum. Hence we will also consider rotating black holes. In the next section, we then 
summarize the relevant features of black holes with Q ^ 0, J ^ 0. 
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(a) (b) 
Figure 7: "White hole" structure in e + e~ annihilation (a) and hadronic collisions (b) 




h h 

Figure 8: Hadronization in hadron-hadron collisions 

5 Charged and Rotating Black Holes 

As mentioned, an outside observer the only characteristics of a black hole are its mass M, 
its electric charge Q, and its spin or angular momentum J. Hence any further observables, 
such as the event horizon or the Hawking temperature, must be expressable in terms of 
these three quantities. 

The event horizon of a black hole is created by the strong gravitational attraction, which 
leads to a diverging Schwarzschild metric at a certain value of the spatial extension R. 
Specifically, the invariant space-time length element ds 2 is at the equator given by 



ds 2 = (1 - 2GM/R) dt 2 - 



1 



1 - 2GM/R 



dr 2 



(64) 



with r and t for flat space and time coordinates; it is seen to diverge at the Schwarzschild 
radius Rs = 2GM. If the black hole has a net electric charge Q, the resulting Coulomb 
repulsion will oppose and hence weaken the gravitational attraction; this will in turn 
modify the event horizon. As a result, the corresponding form (denoted as Reissner- 
Nordstrom metric) becomes 



ds 2 



(1 - 2GM/R + GQ 2 /R 2 ) dt 2 - , _ 2GM/R + GQ2/R2 



dr 2 



(65) 
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For this, the divergence leads to the smaller Reissner-Nordstrom radius 



Rrn = GM (1 + a/1 - Q 2 / GM 2 ) 



(66) 



which reduces to the Schwarzschild radius Rs for Q = 0. The temperature of the Hawking 
radiation now becomes [24, 39] 



T BH (M,Q) =T bh (M,Q) 



4 



Q 2 /GM 2 



1 + Jl- Q 2 /GM 2 ) 2 



(67) 



its functional form is illustrated in Fig. [9l We note that with increasing charge, the 
Coulomb repulsion weakens the gravitational field at the event horizon and hence de- 
creases the temperature of the corresponding quantum excitations. As Q 2 — > GM 2 , the 
gravitational force is fully compensated and there is no more Hawking radiation. 



1.0 



0.5 



T bh (M,Q)/T bh (M,0) 






Q 2 /GM 2 



0.5 1.0 

Figure 9: Radiation temperature for a charged black hole 



In a similar way, the effect of the angular momentum of a rotating black hole can be 
incorporated. It is now the centripetal force which counteracts the gravitational attraction 
and hence reduces its strength. The resulting Kerr metric must take into account that in 
this case the rotational symmetry is reduced to an axial symmetry, and with 9 denoting 
the angle relative to the polar axis 9 = 0, it is (at fixed longitude) given by 

2 / 2GMR \ , R 2 +j 2 cos 2 9 j 2 2 2 2 2 

ds = i 1 - & + j *cos*e ) dt - R 2 -2GMR + j 2 dr ~ {R +J C ° S 9)de ■ (68) 

The angular momentum of the black hole is here specified by the parameter j = J/M; for 
a = 0, we again recover the Schwarzschild case. The general situation is now somewhat 
more complex, since eq. ( 1681) leads to two different divergence points. The solution 



R K = GM (1 4- y/l-j 2 /{GM) 2 ) (69) 

defines the actual event horizon, corresponding to absolute confinement. But the resulting 
black hole is now embedded in a larger ellipsoid 



R E = GM (1 + J I - [j 2 /(GM) 2 } cos 2 9), (70) 
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as illustrated in Fig. [TUJ The two surfaces touch at the poles, and the region between them 
is denoted as the ergosphere. Unlike the black hole proper, communication between the 
ergosphere and the outside world is possible. Any object in the ergosphere will, however, 
suffer from the rotational drag of the rotating black hole and thereby gain momentum. We 
shall return to this shortly; first, however, we note that the temperature of the Hawking 
radiation from a rotating black hole becomes 

T BH (M, J) = T BH (M, 0) V . 

\l + y/l-jy(GM)* 

For a non-rotating black hole, with j = 0, this also reduces to the Hawking temperature 
for the Schwarzschild case. 



ergosphere 










black hole 



Figure 10: Geometry of a rotating black hole 




To illustrate the effect of the ergosphere, imagine radiation from a Schwarzschild black hole 
emitted radially outward from the event horizon. In the case of a Kerr black hole, such an 
emission is possible only along the polar axis; for all other values of 8, the momentum of the 
emitted radiation (even light) will increase due to the rotational drag in the ergosphere. 
This effect ceases only once the radiation leaves the ergosphere. Since the amount of 
drag depends on 9, the momentum of the radiation emitted from a rotating black hole, 
as measured at large distances, will depend on the latitude at which it is emitted and 
increase from pole to equator. 

Finally, for completeness, we note that for black holes with both spin and charge (denoted 
as Kerr-Newman) , the event horizon is given by 



R KN = GM (1 + y/1 - [Q 2 /GM 2 ] - [j'V(GM) 2 ]), (72) 
and the radiation temperature becomes [24, 39] 



4J1- (GQ 2 +f)/(GM) 2 

T BH (M,Q, J) = T BH (M, 0,0) . V - . (73) 

1 + Jl-(GQ*+j*)/(GM)*) 2 + f/(GMY 



The decrease of T BB for Q ^ 0, J ^ expresses the fact that both the Coulomb repulsion 
and the rotational force counteract the gravitational attraction, and if they win, the black 
hole is dissolved. 
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The dependence of a black hole on its basic properties M, Q, J is very similar to the 
dependence of a thermodynamic system on a set of thermodynamic observables. The first 
law of thermodynamics can be written as 

dE = TdS + <f)dQ + oodJ, (74) 

expressing the variation of the energy with entropy S, charge Q and spin J; here denotes 
the electrostatic potential per charge and uj the rotational velocity. The corresponding 
relation in black hole thermodynamics becomes 

dM = T BH dS BH + ®dQ + VldJ, (75) 

where the entropy Sbh is defined as the area of the event horizon, 

obh = ^ • (76) 

The temperature is given by eq. fl75|) . and 

*=^», n =£L (77) 

specify the electrostatic potential $ and the rotational velocity fl. 

The considerations of this section were for spherical black holes. As seen above, such 
objects are in fact equivalent to uniformly accelerating systems. An application to actual 
high energy collisions involves a further assumption. Thermal Hawking-Unruh radiation 
arises already from a single QQ system, as seen above in the discussion of e + e~ anni- 
hilation. If we treat the systems produced in heavy ion collisions as black holes of an 
overall baryon number and/or an overall spin, we are assuming that the collision leads to 
a large-scale collective system, in which each accelerating parton is affected by totality of 
the other accelerating partons. This assumption clearly goes beyond our event horizon 
conjecture and, in particular, it need not be correct in order to obtain thermal hadron 
production. 



6 Baryon Density and Angular Momentum 
6.1 Vacuum Pressure and Baryon Repulsion 

We now want to consider the extension of charged black hole physics to color confinement 
in the case of collective systems with a net baryon number. In eq. ( |67j) we had seen that 
the reduction of the gravitational attraction by Coulomb repulsion in a charged black hole 
modifies the event horizon and hence in turn also the temperature of Hawking radiation. 
The crucial quantity here is the ratio Q 2 /GM 2 of the repulsive overall Coulomb force, 
Q 2 /R 2 , to the attractive overall gravitational force, GM 2 /R 2 , at the horizon. 

In QCD, we have a "white" hole containing colored quarks, confined by chromodynamic 
forces or, equivalently, by the pressure of the physical vacuum. If the system has a non- 
vanishing overall baryon number, the baryon-number dependent interaction will also affect 
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the forces at the event horizon. The simplest instance of such a force the repulsion between 
quarks due to Fermi statistics, but more generally, there will be repulsive effects of the 
type present in cold dense baryonic matter, such as neutron stars. The resulting pressure 
will modify the confinement horizon and hence lead to a corresponding modification of 
the Hawking-Unruh temperature of hadronization. 

By using the conjectured correspondence between black hole thermodynamics and the 
thermodynamics of confined color charges, we translate black hole mass, charge and grav- 
itational constant into white hole energy, net baryon number and string tension, 

{M, Q, G} <-> {E, B, l/2a}. (78) 

Hence eq. (1571) leads us to the relation 

f 4 Jl-2aB 2 /E 2 ) 

T Q (B) = Tq(B = 0) V , ; (79) 

\(1 + s Jl-2aB 2 /E 2 ) 2 J 

for the dependence of the hadronization temperature on the ratio of net baryon number 
B and energy E, with Tq(B = 0) given by eq. (1551) . Its functional form is the same as 
that illustrated in Fig. [91 

It would be interesting to test the prediction (1791) against experimental data; one could 
identify B with the net baryon number per unit rapidity dNs/dy and E with the total 
transverse energy per unit rapidity dET/dy. The reduction of the hadronization temper- 
ature with baryon number could thus occur in two ways. A sufficient decrease of the 
collision energy, e.g. from peak SPS to AGS energy, will strongly reduce dET/dy, while 
dNs/dy is not affected as much. This leads to the known decrease of T(hb) with increas- 
ing n B [44], and it will be interesting to see if the form (1791) agrees with the observed 
behaviour. A second, novel possibility would be to consider hadrochemistry as a func- 
tion of rapidity. At peak SPS energy, dNs/dy remains essentially constant out to about 
y = 2, while dET/dy drops by more than a factor of two from y — to y — 2 [45]. A 
similar behaviour occurs at still lower collision energies. Hence it would seem worthwhile 
to check if an abundance analysis at large y indeed shows the expected decrease of the 
hadronization temperature. 

6.2 Angular Momentum and Non-Central Collisions 

The dependence of Hawking radiation on the angular momentum of the emitting sys- 
tem introduces another interesting aspect for the "white hole evaporation" we have been 
considering. Consider a nucleus-nucleus collision at non-zero impact parameter b. If the 
interaction is of collective nature, the resulting interaction system may have some angular 
momentum orthogonal to the reaction plane (see Fig. 1TT1) . In central collisions, this will 
not be the case, nor for extremely peripheral ones, where one expects essentially just 
individual nucleon-nucleon collisions without any collective effects. 

If it possible to consider a kinematic region in which the interacting system does have an 
overall spin, then the resulting Hawking radiation temperature should be correspondingly 
reduced, as seen in eq. (ITT]) . The effect is not so easily quantified, but simply a reduction 
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Figure 11: Rotating interaction region in non-central AA collision 



of the hadronization temperature for non-central collisions would quite indicative. Such a 
reduction could appear only in the temperature determined by the relative abundances, 
since, as we shall see shortly, the transverse momentum spectra should show modifications 
due to the role of the ergosphere. 

We next turn to the momentum spectrum of the Hawking radiation emitted from a ro- 
tating white hole. As discussed in section 4, such radiation will exhibit an azimuthal 
asymmetry due to the presence of the ergosphere, which by its rotation will affect the 
momentum spectrum of any passing object. At the event horizon, the momentum of all 
radiation is determined by the corresponding Hawking temperature (1711) ; but the passage 
of the ergosphere adds rotational motion to the emerging radiation and hence increases 
its momentum. As a result, only radiation emitted directly along the polar axis will have 
momenta as specified by the Hawking temperature; with increasing latitude 9 (see Fig. 
[121) a. the rotation will increase the radiation momentum up to a maximum value in the 
equatorial plane. 



J 

J i 




(a) (b) 

Figure 12: Transverse plane view of a non-central AA collision 

Hawking radiation from a rotating source thus leads for nuclear collisions quite naturally 
to what in hydrodynamic studies is denoted as elliptic flow. It is interesting to note 
that both scenarios involve collective effects: while in hydrodynamics, it is assumed that 
non-central collisions lead to an azimuthally anisotropic pressure gradient, we have here 
assumed that such collisions lead to an overall angular momentum of the emitting system. 

Concluding this section we emphasize that the results obtained here for the Hawking 
temperatures of systems with finite baryon density or with an effective overall spin depend 
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crucially on the assumption of collectivity. If the different nucleon-nucleon interactions 
in a heavy ion collision do not result in sufficiently collective behavior, the corresponding 
modifications of T q do not apply. In the case of black holes with spin, we moreover have 
no way to relate in a quantitative way centrality and overall spin. Both cases do show, 
however, that such extensions lead to qualitatively reasonable modifications. 



7 Temperature and Acceleration Limits 

We had seen that the underlying confinement dynamics of high energy hadron collisions 
and e + e~ annihilation led to a limit on the acceleration (or the corresponding deceleration) 
in the self-similar hadronization cascade - a limit which can be specified in terms of the 
string tension. In turn, this led to a limiting Unruh hadronization temperature 




Tq ^ J-z-. (80) 



We emphasize that a Hawking-Unruh temperature as such can a priori have any value; 
it is the universal limit on the acceleration that leads to a universal temperature for the 
emitted hadron radiation. 

In the study of strongly interacting matter, temperature limits are well-known and arise 
for an ideal gas of different composite constituents ( "resonances" or "fireballs" of varying 
mass M), if the composition law provides a sufficiently fast increase of the degeneracy 
p(M) with M. If the number of states of a constituent of mass M grows exponentially, 

p(M) ~ M' a exp{bM}, (81) 

with constants a and b, then the grand canonical partition function for an ideal gas in a 
volume V 



N iv • 



V 



(2vr) ; 



-i N 

J dMp(M) J d 3 p exp{-^Jp 2 + M 2 /T] 



(82) 



diverges for 

T>T H = 1/6, (83) 

so that the Hagedorn temperature Th [18] constitutes an upper limit for the temperature 
of hadronic matter. 

In the dual resonance model [46] , the resonance composition pattern is governed by linearly 
rising Regge trajectories, 

a'Ml = n + a 0) n = l,2, (84) 

in terms of the universal Regge slope a' ~ 1 GeV -2 and a constant (of order unity) 
specifying the family (ir, p, ...). For an ideal resonance gas in D — 1 space and one time 
dimension, one then obtains [47] 



1 / Da' 

a = -(£> + !), 6 = 2W— , (85) 
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leading to the temperature limit 




arVo-D- (86) 

In string theory, the Regge resonance pattern is replaced by string excitation modes, 
retaining the same underlying partition structure, with a' = l/2n<j relating Regge slope 
and string tension. Hence we get 

for the corresponding limiting temperature. For D — 4, this coincides with eq. ([IT) 
from [16] and agrees within 20 % with the Unruh temperature flHUl) determined by the 
lowest string excitation alone (cq = 1 in eq. fl58|) ). 

Prior to the dual resonance model, Hagedorn had determined the level density p(M) of 
fireballs composed of fireballs, requiring the same composition pattern at each level [18]. 
The resulting bootstrap condition leads to [48] 



r 



— (21n2-l) 
4 



-1/3 

r , (88) 



where r measures the range of the strong interaction. With r ~ 1 fm, we thus get 
Th — 0.2 GeV for the limiting temperature of hadronic matter. If we identify r with the 
pair production separation x q obtained in eq. ( |54l) . we get 

T H = -^J^ (89) 

X q V Z7T 

and hence again agreement with the hadronic Unruh temperature (IBUI) . 

Hadronic matter as an ideal gas of constituents with self-similar composition spectra 
("resonances of resonances" or "fireballs of fireballs") thus leads to an upper limit of 
the temperature, because the level density of such constituents increases exponentially. 
What does this have to do with the limiting acceleration found in the qq cascade of e + e~ 
annihilation? 

To address this problem, it is useful to recall the underlying reason for the exponential 
increase of the level density in the dual resonance model and the bootstrap model. The 
common origin in both cases is a classical partition problem, which in its simplest form 
[40] asks: how many ways p(M) are there to partition a given integer M into ordered 
combinations of integers? As example, we have for M = 4 the partitions 4, 3+1, 1+3, 
2+2, 2+1+1, 1+2+1, 1+1+2, 1+1 + 1+1; thus here p(M = 4) = 8 = 2 M ~\ It can be 
shown that this is generally valid, so that 

p{M) = ^exp{Mln2}. (90) 

For a "gas of integers", T = l/ln2 would thus become the limiting temperature; the 
crucial feature in thermodynamics is the exponential increase in the level density due to 
the equal a priori weights given to all possible partitions. 
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Returning now to the quark cascade in e + e~ annihilation, we note that the form we have 
discussed above is a particular limiting case. We assumed that the color field of the 
separating qq excites in the first step one new pair from the vacuum; in principle, though 
with much smaller probability, it can also excite two or more. The same is true at the 
next step, when the tunnelling produces one further pair: here also, there can be two 
or more. Thus the e + e~ cascade indeed provides a partition problem of the same kind. 
What remains to be shown are the two specific features of our case: that the dominant 
decay chain is one where in each step one hadron is produced, which provides the constant 
deceleration of the primary quark and antiquark. 

The statistical bootstrap model as well as the dual resonance model lead to self-similar 
decay cascades, starting from a massive fireball (or resonance), which decays into further 
fireballs, and so on, until at the end one has light hadrons. In Fig. [T3~a we illustrate such 
a cascade for the case where the average number k of constituents per step in the decay 
(or composition) partition pattern 

M -> Mn + M 12 + ... + M lk ; M u -> M 21 + M 22 + ... + M 2k ; ... (91) 

is k = 3. In the statistical bootstrap model, k can be can be determined [49]; it is found 
that the crucial feature here is the power term multiplying the exponential increase in eq. 
(IHTT) . For a < 5/2, the distribution in k is given by 

/ s (ln2) fc - x , . 

so that the average becomes 

k = l + 21n2~2.4. (93) 

The dominant decay (~ 70%) is thus into two constituents, with 24 % three-body and 6 
% four-body decays. While in general the fireball mass M could decrease in each step by 
M/k, i.e., by an amount depending on M, the case a < 5/2 is found to be dominated by 
one heavy and one soft light hadron, 

M^Mi + hi; Mx^M 2 + h 2 ; ... (94) 

where hi denotes final hadrons; the pattern is shown in Fig. [13b. Moreover, the three- 
and four-body decays also lead to one heavy state plus soft light hadrons. The decay thus 
provides a uniform decrease of the fireball mass by the average hadron mass or transverse 
energy. 




(a) (b) 
Figure 13: Fireball decay patterns 
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We therefore conclude that the hadronization pattern we had obtained for e + e~ anni- 
hilation is indeed also connected to the same partition problem as the one leading to 
exponential level densities. 

8 Stochastic vs. Kinetic Thermalization 

In statistical mechanics, a basic topic is the evolution of a system of many degrees of 
freedom from non-equilibrium to equilibrium. Starting from a non-equilibrium initial 
state of low entropy, the system is assumed to evolve as a function of time through 
collisions to a time-independent equlibrium state of maximum entropy. In other words, 
the system loses the information about its initial state through a sequence of collisions 
and thus becomes thermalized. In this sense, thermalization in heavy ion collisions was 
studied as the transition from an initial state of two colliding beams of "parallel" partons 
to a final state in which these partons have locally isostropic distributions. This "kinetic" 
thermalization requires a sufficient density of constituents, sufficiently large interaction 
cross sections, and a certain amount of time. 

From such a point of view, the observation of thermal hadron production in high energy 
collisions, in particular in e + e" and pp interactions, is a puzzle: how could these systems 
ever "have reached" thermalization? Already Hagedorn [50] had therefore concluded 
that the emitted hadrons were "born in equilibrium". Given an exponentially increasing 
resonance mass spectrum, it remained unclear why collisions should result in a thermal 
system. 

Hawking radiation provides a stochastic rather than kinetic approach to equilibrium, with 
a randomization essentially provided by the quantum physics of the Einstein-Podolsky- 
Rosen effect. The barrier to information transfer due the event horizon requires that 
the resulting radiation states excited from the vacuum are distributed according to max- 
imum entropy, with a temperature determined by the strength of the "confining" field. 
The ensemble of all produced hadrons, averaged over all events, then leads to the same 
equilibrium distribution as obtained in hadronic matter by kinetic equilibration. In the 
case of a very high energy collision with a high average multiplicity already one event can 
provide such equilibrium; because of the interruption of information transfer at each of 
the successive quantum color horizons, there is no phase relation between two successive 
production steps in a given event. The destruction of memory, which in kinetic equili- 
bration is achieved through sufficiently many successive collisions, is here automatically 
provided by the tunnelling process. 

So the thermal hadronic final state in high energy collisions is not reached through a 
kinetic process; it is rather provided by successively throwing dice. 

9 Conclusions 

We have shown that quantum tunnelling through the color confinement horizon leads to 
thermal hadron production in the form of Hawking-Unruh radiation. In particular, this 
implies: 
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• The radiation temperature Tq is determined by the transverse extension of the color 
flux tube, giving 




in terms of the string tension a. 

• The multiplicity u(s) of the produced hadrons is approximately given by the increase 
of the flux tube thickness with string length, leading to 

i/(s) ~ In y/s, (96) 

where y/s denotes the cms collision energy Parton evolution and gluon saturation 
will, however, increase this, as will early hard production. The universality of the 
resulting abundances is, however, not affected. 

• The temperature of Hawking radiation can in general depend on the charge and 
the angular momentum of the emitting system. The former here provides a baryon- 
number dependence of the hadronization temperature and predicts a decrease of Tq 
for sufficiently high baryon density. The latter provides the basis for the possibility 
of elliptic flow and of a dependence of Tq on the centrality of AA collisions. 

• The limiting temperature obtained in the statistical bootstrap and the dual reso- 
nance or string model arises from a self-similar composition pattern leading to an 
exponentially growing level density. We find that the underlying partition problem 
also leads to the cascade form obtained for hadron emission in high energy collisions, 
so that the dynamic and the thermodynamic limits have the same origin. 

• In statistical QCD, thermal equilibrium is reached kinetically from an initial non- 
equilibrium state, with memory destruction through successive interactions of the 
constituents. In high energy collisions, tunnelling prohibits information transfer and 
hence leads to stochastic production, so that we have a thermal distribution from 
the outset. 

We close with a general comment. In astrophysics, Hawking- Unruh radiation has so far 
never been observed. The thermal hadron spectra in high energy collisions may thus 
indeed be the first experimental instance of such radiation, though in strong interaction 
instead of gravitation. 
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